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Abstract. An acyclic mapping from an n element set into itself is a 
mapping tp such that if tp k (x) — x for some k and x, then (p(x) = x. 
Equivalently, ip — ip l+1 = . . . for i sufficiently large. We investigate the 
behavior as n — > oo of a Markov chain on the collection of such map- 
pings. At each step of the chain, a point in the n element set is chosen 
uniformly at random and the current mapping is modified by replacing 
the current image of that point by a new one chosen independently and 
uniformly at random, conditional on the resulting mapping being again 
acyclic. We can represent an acyclic mapping as a directed graph (such 
a graph will be a collection of rooted trees) and think of these directed 
graphs as metric spaces with some extra structure. Heuristic calcula- 
tions indicate that the metric space valued process associated with the 
Markov chain should, after an appropriate time and "space" rescaling, 
converge as n — > oo to a real tree (K-tree) valued Markov process that 
is reversible with respect to a measure induced naturally by the stan- 
dard reflected Brownian bridge. The limit process, which we construct 
using Dirichlet form methods, is a Hunt process with respect to a suit- 
able Gromov-Hausdorff-like metric. This process is similar to one that 
appears in earlier work by Evans and Winter as the limit of chains in- 
volving the subtree prune and regraft tree (SPR) rearrangements from 
phylogenetics. 



1. Introduction 

A mapping (p from the set [n] := {1, 2, . . . , n} into itself may be repre- 
sented as a directed graph with vertex set [n] and directed edges of the form 
(£, <p(i)), i £ [n]. The resulting directed graph has the feature that every 
vertex has out-degree 1 (with self-loops - corresponding to fixed points - al- 
lowed), and any such graph corresponds to a unique mapping. For example, 
the mapping ip : [18] — ► [18] in Table CD corresponds to the directed graph in 
Figure CD 
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Table 1. A mapping from [18] into itself. 



6 13 17 16 




CD 



Figure 1 . The directed graph corresponding to the mapping 
in Table [TJ 

The directed graph may be decomposed into a number of connected com- 
ponents. Each of these components consists of a single directed cycle (pos- 
sibly a self-loop) plus trees rooted at each vertex on the directed cycle (such 
a tree may be a trivial tree consisting of only the root, meaning that the 
only pre-image of that point is its predecessor on the directed cycle). We 
call such rooted trees the tree components of the graph. 

Aldous and Pitman [AP94] describe a procedure for associating a mapping 
of [n] into itself with a lattice reflected bridge path of length 2n, that is, with 
a function b : {0, 1, . . . , 2n} -> {0, 1, 2, . . .} such that 6(0) = b(2n) = and 
\b(k + 1) — b(k)\ = 1 for < k < 2n. The exact details of the procedure 
aren't important for us. However, we note that a tree component with 
t vertices corresponds to a lattice positive excursion path from with 2t 
steps. Such a segment of path records the distance from the root plus 1 in a 
depth-first-search of the tree component. For example, the tree component 
of size 5 consisting of the vertices {1, 10, 13, 16, 17} in Figure Q] corresponds 
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to the excursion shown in Figure [2] after a suitable translation of the time 
axis. In particular, a tree component that consists of just one point (which 
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Figure 2. The excursion corresponding to the tree compo- 
nent rooted at vertex 10 in Figure [2] with the start of the 
excursion shifted to time 0. 

is necessarily a point on a directed cycle) corresponds to an excursion of the 
form b(k - 1) = 0, b(k) = 1, and b(k + 1) = 0. 

Of course, the mapping cannot be recovered from just the lattice reflected 
bridge path. For one thing, some extra marking of distinguished points of 
the zero set of the lattice path is required to split the lattice path up into 
sub-paths corresponding to components of the directed graph. Once this is 
done, the mapping is uniquely specified by the lattice path up to a relabeling 
of the vertices: that is, if two mappings ip and V correspond to the same 
lattice reflected bridge path, then ip = it o <p> o 7r _1 for some permutation it 
of [n]. Conversely, if ip = it o ip o 7r _1 for some permutation 7r of [n], then 
the lattice path corresponding to ip may be obtained from the lattice path 
corresponding to ip by composition with a bijective map of {0, 1, . . . , 2n} 
that preserves lengths of excursions above all levels. That is, if the lattice 
path corresponding to p has k excursions above some level h, then the same 
is true of the lattice path corresponding to ip. 

Suppose now that a mapping of [n] into itself is chosen uniformly at ran- 
dom from the n n possibilities. This is equivalent to choosing the image of 
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each point of [n] independently and uniformly at random from [n]. The cor- 
responding lattice reflected bridge path is not uniformly distributed. How- 
ever, it is shown in |AP94] that if the lattice reflected bridge path is turned 
into a continuous time process by holding it constant between integer time 
points, time is rescaled by 2n, and space is rescaled by na to produce a 
function from [0, 1] into M+, then this stochastic process with cadlag sam- 
ple paths converges in distribution to twice a standard reflected Brownian 
bridge (that is, twice the Brownian bridge reflected at that goes from posi- 
tion at time to position again at time 1). In particular, the proportion 
of vertices that lie on directed cycles converges to the proportion of time the 
standard reflected Brownian bridge spends at 0, which is, of course, 0, so 
that asymptotically almost all vertices are not roots of tree components. The 
asymptotics of the cyclic vertices jointly with the tree vertices are described 
in |AP 94] using the local time at of the reflected bridge and that paper also 
describes an auxiliary "marking" procedure for describing the joint asymp- 
totics of the the component sizes. Some later results in this same vein may 
be found in [AMP051 \KPQ2~[ IDS971 IDG99L iGLOOl IDG041 IFItn2| . 

A mapping p from [n] into itself is acyclic if the only directed cycles in the 
corresponding directed graph are self-loops. That is, each x £ [n] is either a 
fixed point of ip (so that x is a vertex on a self-loop) or ip k (x) x for any k. 
Equivalently, tp = ip e+1 = . . . for £ sufficiently large. For such a mapping, 
each graph component consists of single tree component with a self-loop 
attached to the root, and no auxiliary marking procedure is necessary to 
recover the mapping up to a permutation from the corresponding lattice 
reflected bridge path. It is not hard to show that if we turn the lattice 
reflected bridge path for a uniformly chosen acyclic random mapping into 
a continuous time process indexed by [0, 1] as above, then the resulting 
process also converges to twice a standard reflected Brownian bridge - as one 
would expect from the observation that the cyclic vertices are asymptotically 
negligible for a uniformly chosen random mapping, 

In this paper we are interested in the asymptotic behavior as n — ► oo of a 
simple Markov chain that randomly evolves an acyclic mapping from [n] into 
itself. At each step of the chain, a point of [n] is chosen uniformly at random 
and the image of this point is re-set to a new image chosen independently 
and uniformly at random from [n], conditional on the resulting mapping 
being acyclic. It is clear for each n that this chain is reversible with respect 
to the uniform distribution on the set of acyclic mappings from [n] into itself 
and that the chain converges to this distribution at large times. 

In terms of the corresponding directed graphs, the chain evolves as fol- 
lows. A directed edge is first chosen uniformly at random and deleted. The 
deleted edge is then replaced by another directed edge with the same ini- 
tial vertex but a uniformly chosen final vertex, conditional on the resulting 
graph having no cycles other than self-loops. Note that the effect of such a 
step is the following. 
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• If the deleted edge is a self-loop, its deletion turns the graph com- 
ponent that contained the edge into a rooted subtree. Otherwise, 
the deletion of the directed edge splits the graph component that 
contained it into two pieces, one of which contains a self-loop and 
the other of which is a rooted subtree. 

• In either case, the addition of the new directed edge either attaches 
the root of the subtree to itself by a self-loop, producing an extra 
graph component, or the new directed edge attaches the root to a 
vertex chosen uniformly outside the subtree (possibly to a vertex 
outside the subtree but within the same former graph component). 
All such possibilities are equally likely. 

The effect on the corresponding lattice bridge path is to remove an excur- 
sion above some level, insert a suitable time-space translation of it at some 
time point in the lattice bridge path outside the excursion, and then close 
up the gap left by the removal (more precisely, this transformation may need 
to be followed by a bijective map of {0, 1, . . . , 2n} that preserves lengths of 
excursions above all levels because of the way that the labeling of vertices 
in the directed graph is used to construct the corresponding lattice bridge 
path). 

In order to understand the asymptotic behavior of this sequence of chains 
as n — > oo, we need to embed the state space of each chain into a common 
state space that will also be the state space of the limit process. 

To begin with, we erase all of the self-loops in the directed graph corre- 
sponding to an acyclic mapping. This produces a forest of subtrees rooted 
at vertices that were formerly on self- loops. We then connect the roots of 
these subtrees by directed edges to a single adjoined point to produce a tree 
rooted at the adjoined point. Keeping in mind the rescaling identified by 
Aldous and Pitman, we think of this tree as a one-dimensional cell complex 
by regarding each edge as a segment of length n~2 . We thus have a metric 
space with a distinguished base point (the root). This pointed metric space 
is an instance of a rooted compact real tree (R-tree): see Section [2] for the 
precise definition of a IR-tree - for the moment, all that is important for 
explaining our results is that a M-tree is a metric space that is, in some 
sense, "tree-like". We regard two rooted compact M-trees as being equal if 
one can be mapped into the other by an isometry that preserves the root. If 
two mappings <p and ij) are related by a relabeling ifi = ir o ip o 7r _1 for some 
permutation ir of [n], then they correspond to the same rooted compact 
R-tree. 

Before we continue with the motivation of our results, we need to indicate 
how the rooted compact M-tree associated with a mapping from [n] into itself 
may be constructed directly from the corresponding lattice reflected bridge 
path. We begin by introducing some general notation that is useful later. 
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Definition 1.1. Write C(R_|_,R+) for the space of continuous functions from 
R+ into R+. For / £ C(R+,R+), put 

C(/) := inf{s > : f(t) = for all t > s} 

with the usual convention that inf = oo. The set of positive bridge paths 
is the set Q + C C(R + ,R + ) given by 

/(0) = 0, < C(/) < oo, 



n+ - Y G C '( R +' R +) : f(t) > o for < t < C(/). 
For £ > 0, set Q e + : = {/ £ 0+ : £(/) = £}. 

We associate each / £ fii. with a compact metric space as follows. Define 
an equivalence relation ~j on [0, 1] by letting 

m ~/ u 2 , iff /(tii) = inf /(it) = /(u 2 ). 

Consider the pseudo-metric djy on [0, 1] defined by 

d T Au 1 ,u 2 ) := f(ui) - 2 inf /(-u) + /(u 2 ). 

This pseudo-metric becomes a true metric on the quotient space 1/ := 
[0,1]/^^. The resulting metric space is compact and is an instance of a 
rooted compact R-tree if we define the root to be the image of under the 
quotient map. 

Suppose that the function / £ is obtained by first linearly interpolat- 
ing the lattice reflected bridge path associated with an acyclic mapping ip 
of [n] into itself to produce a function in Q^_ n and then rescaling time by 2n 

and space by wa. The corresponding pointed metric space Tf is the rooted 
compact R-tree associated with tp that we described above. 

Any metric space of the form Tf for / £ f2+ has two natural Borel measure 
on it. Firstly, there is the "uniform" probability measure it, given by the 
push-forward of Lebesgue measure on [0, 1] by the quotient map. We call 
this measure the weight on Tf . Secondly, there is the natural length measure 
HTf, which is the one-dimensional Hausdorff measure associated with the 
metric dr, restricted to points of Tf that are not "leaves" (see Section [2] for 
a more precise definition). When / is associated with a map of [n] into itself 
as above, then fiT f is just the "Lebesgue measure" on the cell complex Tf 

that assigns mass to each edge of Tf (recall that we have rescaled so 
that each edge has length n~2). 

Now, if we speed up time by a factor of nz in our Markov chain for 
evolving mappings of [n] into itself and look at the corresponding rooted 
compact R-tree- valued process, then it is reasonable at the heuristic level 
that we should obtain in the limit as n — * oo a continuous time Markov 
process with the following informal description. The state space of the limit 
process is the space consisting of rooted compact R-trees T equipped with 
a probability measure vt- we call such objects weighted rooted compact 
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M-trees. We note that, as in the special case of M-trees of the form Tf 
for / G an arbitrary compact M-tree has a canonical length measure /it 
given by the restriction of the one-dimensional Hausdorff measure associated 
with the metric to the set of points that aren't leaves. The process evolves 
away from its state at time by choosing a point (t, v) at rate dt (g) fix(dv) 
in time and on the current tree T, and at time t the subtree above v (that 
is, the subtree of points on the other side of v from the root) is re-attached 
at a point w chosen according to vt (conditional on w being outside the 
subtree). 

In general, the measure [it may have infinite total mass. For example, if 
/ G Oi. is chosen according to the distribution of standard reflected Brown- 
ian bridge, so that is the rooted compact M-tree that arises from a limit 
as n —* 00 of uniform acyclic random mappings of [n] into itself, then fj,x 2 * 
almost surely has infinite total mass. Consequently, the above specification 
of the dynamics of the limit process does not make rigorous sense for general 
weighted rooted compact M-trees. The aim of this paper is to use Dirich- 
let form methods to construct a suitably well-behaved Markov process with 
evolution dynamics that conform to the heuristic description. 

We do not, however, obtain a convergence result. The limit process has no 
obvious Feller-like properties and it is not clear how to define its dynamics for 
all starting points (as opposed to almost all starting points with respect to 
the symmetrizing measure, which is all the Dirichlet form approach provides) 
in such a way that, say, martingale problem methods might be used to 
establish convergence. 

The process we construct is somewhat similar to the process constructed 
in |EW06] as a limit a natural chain based on the subtree prune and re- 
graft (SPR) tree rearrangement transformations from phylogenetics. Both 
processes involve the relocation of a subtree whose root is chosen according 
to the length measure on the current tree. However, the state space of the 
process in [EW06J consists of weighted unrooted compact M-trees, whereas 
we work with weighted rooted compact M-trees and the root plays a cru- 
cial role in defining the dynamics. The symmetrizing measures are, as a 
consequence, rather different: the measure in [EW06J is the distribution of 
the Brownian continuum random tree, which is the M-tree "inside" twice a 
standard Brownian excursion, whereas our symmetrizing measure is the dis- 
tribution of the M-tree "inside" twice a standard reflected Brownian bridge. 
However, many of the steps in the construction are quite similar so we omit 
several arguments and simply refer to the analogous ones in [EW06J. 

We note that Markov processes with reflected bridge paths as their state 
space and continuous sample paths have been studied in [Zam03l IZam02t 
IZamOl] . These processes are reversible with respect to the distribution of a 
Bessel bridge of some index. 
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2. Weighted IR-trees 

Definition 2.1. A metric space (T,d) is a real tree (K-tree) if it satisfies 

the following axioms. 

Axiom 0: The space (T, d) is complete. 

Axiom 1: For all x,y £ T there exists a unique isometric embedding 
(j> X) y : [0,d(x,y)] -> T such that (j> X) y(0) = x and (j> x>y (d(x,y)) = y. 
Axiom 2: For every injective continuous map if) : [0, 1] — > T one has 
^([0,1])= ^(o)^(i) ([0,d(^(0),^(l))]). 

Axiom 1 says simply that there is a unique "unit speed" path between 
any two points x and y. We write [x, y] for the image of this path and call 
it the segment with endpoints x and y. Axiom 2 implies that the image of 
any injective path connecting two points x and y coincides with the segment 
[x, y], and so such a path may be re-parameterized to become the unit speed 
path. Thus, while Axiom 1 is satisfied by many other spaces such as IR rf 
with the usual metric, Axiom 2 captures the essence of "treeness" and is 
only s atisfied by R d when d = 1. See |Dre841 IDT961 IDMT951 IDMT961 
ITer97l IChiOlj for background on IR-trees. In particular, [ChiOl] shows that 
a number of other definitions are equivalent to the one above. See also 
[Eva06| , where much of this material is synthesized and combined with other 
material on probability on IR-trees. 

We define the r]-trimming, Rn(T) of a compact IR-tree (T,d) for n > 
to be the set of points x G T such that x belongs to a segment [y, z] with 
d(x, y) = d(x, z) = rj — see Figure El The skeleton of (T, d) is the set 
T° := {J v> q Rr)(T). Thus x G T° if x &]y, z[ for some y, z. The leaf set of 
(T, d) is the set T\T°. The length measure on T is the cr-finite measure /it 
on the Borel cr-field B(T) given by the trace onto T° of the one-dimensional 
Hausdorff measure associated with d. Equivalently, /jlt is the unique measure 
concentrated on T° such that ht([x, y\) = d(x, y) for all x,y £ T (see Section 
2.4 of [EPW06] or Section 2 of |EW06] ). 

In the following, we are interested in compact IR-trees (T, d) equipped 
with a distinguished base point p € T (called the root) and a probability 
measure v on the Borel cr-field B(T) (called the weight). We call such objects 
weighted rooted compact ~R-trees. We say that two weighted rooted compact 
IR-trees ( X, dx , px , vx ) and (Y, dy , Py, vy) are weighted rooted isometric if 
there exists a bijective isometry $ between the metric spaces (X, dx) and 
(Y,dy) such that &(px) = Py and the push-forward of vx by $ is vy> that 
is, 

vy = := vx o 

The property of being weighted rooted isometric is an equivalence relation. 
We write T wr for the collection of equivalence classes of weighted rooted 
compact IR-trees. 

In order to define a metric on T wr , we first recall the definition of the Pro- 
horov distance between two probability measures (see, for example, [EK86] ) . 
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Figure 3. A R-tree T and its jy-trimming R V (T). The R-tree 
T consists of both the solid and dashed segments, whereas 
the R-tree RqiT) consists of just the solid segments. 

Given two probability measures a and (5 on a metric space (X, d) with the 
corresponding collection of closed sets denoted by C, the Prohorov distance 
between them is 

d P (a,/3) := mi{e > : a{C) < (3(C £ ) + e for all C G C}, 

where C £ := {x G X : ini ye c d(x,y) < e}. The Prohorov distance is a 
metric on the collection of probability measures on X. 

We are now in a position to define the weighted rooted Gromov-Hausdorff 
distance between the two weighted rooted compact R-trees (X, dx, px, vx) 
and (Y,dy,pK,&y). 

For e > 0, let F XY denote the set of Borel maps / : X — > Y such that 
f{px) = PY and 

sup{\d x (x',x") - d Y (f(x')J(x"))\ : x',x" e X} < e, 
and define Fy x similarly. Put 
A GH w r (X,Y) 

exist / G F XY ,g £ F YX such that 1 
dp(f*vx,VY) < £, dp(v x ,g*VY) < e f ' 
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Note that the set on the right hand side is non-empty because X and Y 
are compact, and hence bounded in their respective metrics. Note also that 
Agh wi (X, Y) only depends on the weighted rooted isometry classes of X 
and Y. 

It turns out that the function Aqh wi satisfies all the properties of a metric 
except the triangle inequality. To rectify this, put 

dGR*x(X,Y) := inf j Y] A G H" (Zj , Z i+X ) i j , 

where the infimum is taken over all finite sequences of weighted rooted com- 
pact R-trees Z\, . . . Z n with Z\ = X and Z n = Y (the exponent \ is not par- 
ticularly important, any sufficiently small number would suffice). Note again 
that c?GH wr (X, Y) only depends on the weighted rooted isometry classes of 
X and Y. 

From now on, we think of Agh wi and dGH wr as being defined on T wr x T wr . 
Parts (i) and (ii) of the following result are analogous to Lemma 2.3 of 
[EW06], part (iv) is analogous to Proposition 2.4 of [EW06j, part (v) is a re- 
statement of Lemma 2.6 of [EW06], and part (vi) is analogous to Theorem 
2.5 of [EW06]. The results in [EW06] are for M-trees with weights but 
without roots, but the addition of roots does not present any new difficulties 
(c/. the passage from M-trees without weights or roots to M-trees without 
weights but with roots in Section 2.3 of [EPW06]). The space T in part (iv) 
is the collection of isometry classes of compact M-trees (without weights or 
roots) and we refer the reader to Section 2.1 of [EPW06J for the definition 
of the associated Gromov-Hausdorff distance c?gh- 

Proposition 2.2. (i) The map Agh wi has the properties: 

(a) A GR ^(X,Y) = if and only ifX = Y, 

(b) A GH wr(X,Y) = A GH wr(y,X). 

(ii) The map c?GH wr is a metric on T wr . 

(iii) For all X,Y £ T wr , 

iA GH wr(X,y)i < d Gll ™(X,Y) < A GH wr(X,y)3. 

(iv) A subset D of (T wr ,G?GH wr ) is relatively compact if and only if the 
subset E := {(T,d) : (T,d,p,u) £ D} o/(T, (2gh) is relatively com- 
pact. 

(v) A subset E of (T, don) is relatively compact if and only if 

sup{fiT{R n {T)) : T G E} < oo 

for all rj > 0. 

(vi) The metric space (T wr ,dGH wr ) is complete and separable. 

We note that an extensive study of spaces of metric spaces equipped with 
measures is given in |Stu06a | IStu0 6b|. and the theory of weak convergence 
for random variables taking values in such spaces is developed in [GP W06] . 
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3. Trees and continuous paths 

Definition 3.1. The space of positive excursion paths is the set C 
0+ C C(M + ,M + ) given by 

o ,_ If r CCu M ) ■ /(O) = 0, < C(/) < oo, 
.- j/ G G(K +) M+) . > Q for Q < t < c(/) 

For £ > 0, set Q e ++ := {/ G Q++ : C(/) = (}■ 

The following result is a slight generalization of Lemma 3.1 in [EW06]. 
The latter result was for the special case of M-trees constructed from positive 
excursion paths rather than general positive bridge paths. The proof goes 
through unchanged. 

Lemma 3.2. For each f G Q^_, the metric space (Tf,dr.) is a compact 
M-tree. 

We root a M-tree {Tf,dx s ) coming from a positive bridge path in / £ Q\ 
by taking the root to be the point corresponding to £ [0, 1] under the 
quotient map. We equip (Tf,d,T f ) with the weight z/y„ given by the push- 
forward of Lebesgue measure on [0, 1] by the quotient map. 

For a positive bridge path / G O^, we identify the length measure /j,?, 
on the associated compact M-tree (Tf , drp. ) as follows (the discussion is es- 
sentially the same as that in Section 3 of [EW06] which considered M-trees 
coming from positive excursion paths). For a > 0, let 



(3.1) G(f,a):-- 



f(s) = a and, for some t > s, 
s £ [0, 1] : /(r) > a for all r e]s, t[, 
f(t)=a. 



denote the countable set of starting points of excursions of the function / 
above the level a. Then, the length measure fiT f is the push-forward of the 
measure 

St 

teg(f,a) 



/•CO 

(3.2) m f := da V 



by the quotient map, where 5t is the unit point mass at t. 
Alternatively, write 

T(/) :={(«, a) : s G]0, 1[, a G [0, f(s)[} 

for the region between the time axis and the graph of /, and for (s, a) G T(/) 
denote by 

(3.3) s(f,s,a) := sup{r < s : f(r) = a} 
and 

(3.4) s(f,s,a):=mf{t>s:f(t) = a} 
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the start and finish of the excursion of e above level a that straddles time s. 
Then, 

(3.5) m/= / ds<g>da— 1 \ S 8(j,s,a)- 

Mf) s (f,s,a) - s(f,s,a) - KJ ' ' ' 

4. A PATH TRANSFORMATION CONNECTING REFLECTED BROWNIAN 
BRIDGE AND BROWNIAN EXCURSION 

Write P + for the law of the standard Brownian bridge reflected at that 
goes from at time to at time 1. Write P ++ for the law of standard 
Brownian excursion. Thus, P+ is a probability measure on and P++ 
is a probability measure on We show in this section how various 

computations for P + can be reduced to computations for P ++ using a result 
of Bertoin and Pitman. 

Given / G Q+, put 

'limsup £i0 £ / [o t] ds l{/(s) < e}, if limsup ei0 £ J [Q £] ds l{/(s) < 
L(t;f):=l <oo, 

k 0, otherwise, 

for < t < £, and set L(t; f) = L(£; f) for t > i. 

Denote by the subset of Q e + consisting of functions / with the prop- 
erties: 

• the closed set {t G [0,£] : f(t) = 0} is perfect (that is, has no isolated 
points) and has Lebesgue measure zero; 

• for < t < £, 

L(t;/) = Iiml / dsl{f(s)<e}; 
£ i° 2 e J[o,t] 

• the function t ^ L(t; f) is continuous; 

• the set of points of increase of the function t *— ► L(t; f) coincides 
with {t G [0,£] : f(t) = 0}. 

Note that if / G then L(-; /) is not identically (indeed, L(-; /) has 
as a point of increase). Of course, P + ($l 1 |_) = 1. 
For / Gfi' + , set 

[/(/) := sup |o < t < £ : L(i; /) < h(£; /)| 

and put 

f L(t;f), 0<t <U(f), 

K^(t; f) := { L(£; f) - L(t; /), U(f) <t<£, 
0, t>£. 
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For /erf + and u £ [0,£], set 

mint< s <„ f(s), < t < u, 
K*~(t',f,u) ■= I mm u < s < t f(s), u<t<£, 
0, t>£. 
The following result is elementary and we leave the proof to the reader. 
Lemma 4.1. Fix a function f £ fq_. Set 

e = K-*(-;f) + f. 

Then, e £ 

f = K^(-;e,U(f)). 

The next result, which is Lemma 3.3 of [B P94] . says that under P + the 
path- valued random variable / i— > K~*(-;f) + / has law P++, the random 
variable / i— ► U{f) is uniformly distributed on [0, 1], and these two random 
variables are independent. 

Proposition 4.2. For any Borel function F : Q\_ x [0, 1] — > R + , 

/ P + (d/)F(i^(-; /) + /,[/(/)) = / P ++ (de) / duF(e,u). 
J ■/ J[o,i] 

In order to apply Proposition 14.21 we need to understand for a fixed 
positive bridge path / £ Cl^_ how the measure mj of (|3.2|) or (|3.5|) is related to 
the analogous measure for the associated positive excursion path K^(-; f) + 
/ G & ++ . 

Definition 4.3. For e £ and u £ [0, 1], write 

r*(e, u) := {(s, a) £ T(e) : u ^ [s(e, s, a), s(e, s, a)]} 

for the set of points in T(e) such that the corresponding straddling sub- 
excursion does not straddle the time u. 

Lemma 4.4. Fix f £ 0+. Set e = K~*(-\ f) + f £ so that 

r*(e, [/(/)) = {(s, a) : * £]0, 1[, K^(s; f) < a < K^(s; f) + f(s)}. 

Define a bijection £ : T(f) — > r*(e, U(f)) by setting 

£(s,a) := (s,a + ir _+ (s;/)). 

JTie map £ is a measure-preserving bijection between the set T(f) equipped 
with the measure 

ds <g) da — — 

s{f,s,a) -s{f,s,a) 

and the set T*(e, U(f)) equipped with the measure 

ds (g) da 



s(e, s, a) — s(e, s, a) 
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Proof. Decompose the open set {t G [0, 1] : f(t) > 0} into a countable union 
of intervals A k , k G N. Set B k = {(s,a) G T(f) : s G kN, and 

C fc = {(a, a) G r*(e, [/(/)) : s G A k }, k G N. We have A([0, 1] \ \J k A k ) = 0, 
where A is Lebesgue measure. Thus, A ig) A(r(/) \ Ufc^fc) = ana ^ ^ ® 

A(r*(e, [/(/)) \U fc c fc )=o. 

The function t i— ► L(t; f) is constant on each of the sets ^4/%, and so the 
same is true of the function t t— > K~^(t; /). Write for this constant. The 
function £ maps bijectively into C k and the restriction of £ to -Bfc is the 
translation (s, a) i— * (s, a + c^). 

Therefore, £ is a measure-preserving bijection between the set T(f) equipped 
with the measure ds<8>da and the set T*(e, U (/)) equipped with the measure 
ds £g> da. 

It remains to note that if, for some (s, a) G we write a) = (s, a'), 

then we have s(f,s,a) = s(e,s,a') and s(f,s,a) = s(e,s,a'), so that, in 
particular, s(f, s, a) — s(f, s, a) = s(e, s, a') — s(e, s, a'). □ 



Remark 4.5. Assume that / G For a > 0, recall the definition of G(f, a) 
from flHZE]) . For u G [0, 1] and e G 0+ + put 



^*(e,o,u) := 



e(s) = a and, for some t > s, 

r „ e(r) > a for all r Gls,tf, 

s G 0, 1 : , v 

e(t) = a, 

u ^ [a,t]. 



That is, G*(e,a,u) is the countable set of starting points of excursions of e 
above the level a that don't straddle the time u. A consequence of Lemmal4T 
is that the measure mj coincides with the measure 



da 



E 



St. 



t&g*(K^(-J)+f,a,U(f)) 

As explained in the Introduction, the dynamics of the process we wish to 
construct involves "picking" a point v in a rooted compact M-tree (T, dr, px) 
according to the length measure \it and then re-rooting the subtree above 
v (that is, the subtree consisting of points of x G T such that v G [pt,x[) 
at a new location w. When T = Tt for some / G f2+, this re-rooting of a 
subtree corresponds to a rearrangement of / by relocating an excursion of 
/ above some level. We introduce the following notation to describe such 
rearrangements. 

Definition 4.6. For / G n\ and (s,a) G r(/), define / s ' a G 0++ and 
f s > a G f2 + , by 



f s ' a (t) 



/(«(/, s,a) + t) - a, < t < s(f, s, a) - s(f, s, a), 



0, 



t > s{f,s,a) -s(f,s,a), 
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/ s ' a (t) := 



and 

'/(*), 0<t<s(f,s,a), 
J(t + s(f,s,a) -s(f,s,a)), t > s(f,s,a). 

That is, f s ' a is the sub-excursion of / that straddles (s, a) shifted to start 
at position at time 0, and f s ' a is / with the sub-excursion that straddles 
(s, a) excised and the resulting gap closed up. 



U(f,u,s,a) 



Definition 4.7. For / G u G [0, 1], and (s,a) G T*(f,u), put 

j u, < u < s(f, s, a), 

\u - s(f, s, a) + s(f, s, a), s(f, s, a) < u < 1. 

By definition of T* (/,«), the point u belongs to the set 

[0,s(f,s,a)[U]s(f,s,a),l] 

of length 

£(/«•«) = 1 - C(r> a ) = 1 - (*(/, S , a) - s(f, s, a)), 

and U(f,u,s,a) is where « is moved to when we close up the gap to form 
the interval [0,C(/ S ' a )]- 

The following result is immediate from Lemma l4.1l and Lemma 14.41 

Corollary 4.8. Fix f G Set e = K^(-; f) + / G Then, for any 

Borel function F : S7 + x £1+ — > IR +; 

ds ® da ^7 ^T7 V F (^' a ' /' ,a ) 

r(/) s(f,s,a)-s{f,s,a) 

j ds ® da -F(e s ' a , e s ' a , C/(e, [/(/), s, a)). 

Ir*( e ,u{f)) s(e, s, a) - s(e, s, a) 



5. Standard Brownian excursion and length measure 

We first recall a result (Proposition 15.21 below) that appears as Corollary 
5.2 in [EW06]. It says that if we pick an excursion e according to the 
standard excursion distribution IP++ and then pick a point (s,a) G T(e) 
according to the cr-finite measure 

ds (8) da — r -. r 

s(e, s, a) — s{e, s, a) 

so that the time point s(e, s, a) is picked according to the cr-finite measure 
m e , then the following objects are independent: 

(a) the length of the excursion above level a that straddles time s; 

(b) the excursion obtained by taking the excursion above level a that 
straddles time s, turning it (by a shift of axes) into an excursion e s ' a 
above level zero starting at time zero, and then Brownian re-scaling 
e s,a to produce an excursion of unit length; 
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(c) the excursion obtained by taking the excursion e s ' a that comes from 
excising e s ' a and closing up the gap, and then Brownian re-scaling 
e s ' a to produce an excursion of unit length; 

(d) the starting time s(e, s,a) of the excursion above level a that strad- 
dles time s rescaled by the length of e s,a to give a time in the interval 
[0,1]. 

Moreover, 

• the length in (a) is "distributed" according to the cr-finite measure 

1 dr 
2 V ^F v /(i-r> 3= ' rG[ °' 1J; 

• the unit length excursions in (b) and (c) are both distributed as 
standard Brownian excursions (that is, according to 

• the time in (d) is uniformly distributed on the interval [0, 1]. 

Definition 5.1. For c > 0, let S c : — > £1+ be the Brownian re-scaling 
map defined by 

s c f ■= VZfWc). 



Proposition 5.2. For any Borel function F : [0, 1] x f&++ x fL 

r ++ (d») / *° - f(^0 ,**,** 

h(e) s(e, s, a) - s(e, s, a) V C(e s ' a ) 

= / cb-i= / dr [ P ++ (d e ')0P ++ (de")F( W ,5 r e',5 1 _ r e"). 

J[o,i] 2V27T J[ ,i] VI 1 - r ) r J 

With Proposition 14.21 and Corollary 14.81 in mind, we want to obtain an 
analogous result with T(e) replaced by T*(e,u), where u is picked uniformly 
from [0, 1] . 

Corollary 5.3. For any Borel function G : [0, 1] x [0, 1] x fi ++ x Q ++ — > M +; 

ds da 



[ du [ I 


P ++ (de) / 


/[0,1] J 


Jr*{e,u) 



' U(e,u,s,a) s(e,s,a) 
' C(e s < a ) ' C(e s ' a ) 



du / dv-±= [ dr Jl-^ / P ++ (de')®P ++ (de") 
[o,i] 7[ ,i] 2V2^7[o,i] V H 

x G(u,v,S r e' ,Si- r e"). 



Proof. For v,r £ [0, 1] and u e [0, (1 — r)i>[U](l — r)i> + r, 1], put 

< u < (1 — r)v, 



U (u, v, r) 



l-r> 

fEf, (l-r)v + r<«<l. 
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From Proposition 15.21 we have 

ds (g) da 



J P ++ (de) J 



x G 



r*(e,u) s(e,s,a) - s(e,s,a) 
U(e,u, s,a) s(e,s,a) 



((e s > a ) ' C(e s ' a ) 

= /*P ++ (de) / ds ® d( * -l{u(£[s(e,s,a),s(e,s,a)]} 

J Jr(e) s{e,s,a) - s{e,s,a) 

( U(e,u,s,a) s(e,s,a) \ 

= / dv -±= I ^ 3 / P ++ (de') P ++ (de") 

J[o,i] 2V27T J[o,i] y^l - r)H J 

x l{w, ^ [(1 — r)v, (1 — r)u + r]} 
x G ^ (u, v, r), v, S r e , 5i_ r e"^ . 

The change of variable ui = v, r) gives 

/ *,!{. i [(1 - r),, (1 - r)„ + r])G (&(„,„, r),,,^,^ 

i[0,l] v 7 

= (1 — r) / dw G (u>, v , S r e' , Si- r e") , 
7(0,1] 

and the result follows. □ 



6. A SYMMETRIC MEASURE ON X 

Definition 6.1. Fix a function / £ Ol and suppose that v £ (?(/, a) is the 
starting point of an excursion of / above some level a. Write 

5(f,v) :=mf{t>v:f(t) = a} 

for the time at which the excursion finishes. Thus, s(/, s,a) = v and 
s(f,s,a) = S(f,v) for any s £]v,5(f,v)[. Define e^ ,v € Q++ by 

gf,v . = //(< + - < * < V - 5 Ui V), 

\0, t>v-5(f,v). 

That is, e*'" is the result of taking the excursion starting and ending at times 
v and 5(f, v), respectively, and shifting the time and space axes to obtain an 
excursion that starts at position at time 0. Given w € [0, 1] \ [v,5(f,v)], 
denote by f v,w G the path defined as follows. If w > v (so that w > 
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<K/»), then 
f v ' w (t) ■= < 



(f(t), 0<t<v, 

f{t-v + 5{f,v)), v<t<v-6(f,v)+w, 

e f ' v (t -(v- 6(f, v) + w)) + f(w), v - 5(f, v)+w<t<w, 

[f(t), t>w. 



If w < v, then 



f/(t), 0<t<u;, 

e f ' v (t-w) + f(w), w<t<w-v + S(f,v), 

f(t + v- 6(f, v)), w-v + 5(f, v)<t< 6(f, v), 

lf(t), t>5(f,v). 

In other words, the excursion of / starting at time v is first moved so that 
it starts at w and then the resulting gap left between times v and 5(f,v) is 
closed up - see Figure [H 




f v,w 




Figure 4. The transformation taking the path / to the path 
f v,w when w > v. The figure for w < v is similar. 
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Definition 6.2. Define a kernel k + on by 

K+(f,B):= [°°da Yl i ,,,, , v/ d«,l(re5). 

That is, a starting point v of an excursion is chosen according to the measure 
rrif corresponding to length measure /j,T f on the M-tree associated with /, 
this excursion is then relocated so that it starts at a uniformly chosen point 
w G [0, 1] \ [v,6(f,v)], and finally the resulting gap is closed up. Define a 
measure J + on f2+ x by 

J + (d/',d/") :=P + (d/')/€ + (/',d/"). 
Proposition 6.3. The measure JT+ is symmetric. 

Proof. Given e',e" G v G [0,1], and r G]0, 1], define e°(-; e', e", v, r) G 

n\ + by 

e (i;e',e>,r) 

'<Si_ r e"(i), < i < (1 - r>, 

5i_ r e"((l - r)v) + S r e'(t - (1 - r)«), (1 - r)v < t < (1 - r)u + r, 
Si-re"(t- r), (l-r)v + r <t < 1. 

That is, e°(-; e' , e" , v, r) is the excursion that arises from Brownian re-scaling 
e 1 and e" to have lengths r and 1 — r, respectively, and then inserting the 
re-scaled version of e' into the re-scaled version of e" at a position that is a 
fraction u of the total length of the re-scaled version of e". 
Also, for u G [0, 1] set 

(1 — r)u, < u < v, 
r + (1 — r)u, i> < u < 1, 



£/ (u, v, r) :- 



so that U(u,v,r) belongs to the set [0, (l — r)v[U](l — r)v+r, 1] for Lebesgue 
almost all u G [0, 1] and the push-forward of Lebesgue measure on [0, 1] by 
the map u t— > U(u, v, r) is the uniform distribution on this union of two 
intervals. 

Define a measure on [0, 1] x [0, 1] x n\ + x n\ + by 
J++ (du* ,du**, de* , de** )G(u* , u** , e* , e** ) 

:= / du ® dv ® du>— ^= / dr\ \ P++(de') P ++ (de") 

7[o,i]3 2V2¥7 [0 ,i] V r 3 y 

x G ^£7 (u, v, r), U(u, w, r), e°(-; e', e", u,r), e°(-; e' , e" , w, r)J 

for any non-negative Borel function G. 

Clearly, the measure J ++ is preserved by pushing it forward with the 
map (u* , u** , e* , e**) i— » (u** , u* , e** , e*). Also, it follows from Lemma [47Tj 
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Proposition H21 Corollary 14.81 and Corollary 15.31 that the measure JT+ is the 
push-forward of the measure J++ by the map 

[u ,u ,e ,e ) i— ► (K (-;e ,u ),K (-;e ,u )), 

and the result follows. □ 

By construction, the measure J+ is concentrated on pairs (/', /") G fii. X 
Q]_ such that /" is obtained from f by the re-location of an excursion. If 
we shift the starting point of this excursion in space and time to the origin 
to obtain an element of then the u-finite law of this shifted excursion 
is 



r 00 / 

J da 1 [z f ' V G 



veG(f,a) 

ds <g) da 
r(/) 



1 



s{f,s,a) -s(f,s,a) 



1 l f s,a e 



Informally, this is the law of the excursion under J+, but we note while J+ 
is concentrated on pairs (/', /") of the form (/, f v,w ) for some v, w G [0, 1], 
the value of v and the corresponding excursion &> v cannot be uniquely 
reconstructed from (/', /"). Arguing as in the proof of Proposition 16.31 this 
law is given by 



1 



2V2vr y [0 ,i 



dr 



1 — r. 



{e G n\ + : S r (e) G •} . 



We need the following properties of this law. 
Proposition 6.4. (i) For < t < 1, 



J P + (d/) I" da Yl HC^n > t) = -j= ( + arcsin (ft) - \ j , 
veS(f,a) V / 



and hence 



(df) / da Y (C(^)) 5 



veG(f,a) 



1 

7T2 
16 y/2 



(ii) For x > 0, 



//*oo < ^~ > poo / 2 

P+(d/) da V 1 (max(e^) > x) = V / dz exp ( -1- 



and hence 



r + (d/) ^ do £ (ma*(^)) - ^ 



ue5(/,a) 
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Proof, (i) By the remarks prior the statement of the proposition, the quan- 
tity in the first claim is 

2\/2^ I[t,i] ^ V r 3 ' 

and a straightforward integration shows that this has the stated value. The 
second claim follows by an equally straightforward integration by parts, 
(ii) Again by the remarks prior the statement of the proposition, the quantity 
in the first claim is 

—]= I dr A/— I e G n ++ '■ max(e) > 4= 1 • 
2V2W[o,i] V r3 1 V^) 

From Theorem 5.2.10 in [Kn i81j . we have that 

oo 

P++ {e G n\ + : max(e) > y) = 2 J^(4raV - 1) exp(-2nV), 

n=l 

and an integration establishes the claim. 

An integration by parts shows that the quantity in the second claim is 




as required. □ 

7. A DlRICHLET FORM 

Recall that any / 6 £1+ is associated with a R-tree (Tf , d,T f ) that arises 
as a quotient of [0, 1] under an equivalence relation defined by /. Moreover, 
we may equip this M-tree with the root pr f that is the image of under the 
quotient and the weight that is the push-forward of Lebesgue measure 
on [0, 1] by the quotient map. 

Definition 7.1. Define the probability measure P on T wr to be the push- 
forward of the probability measure P+ on tt\ by the map 

/ i ^ (T 2 f,d T2f ,pT 2f ,u T2f ). 

Define the measure J on T wr x T wr to be twice the push-forward of J+ by 
the map 

(/',/") ^ ( (T 2 fi , dr afl , PT 2f , , vt 2S , ) , (T 2f » , dT 2f „ , pr 2f „ , vt 2S „ ) ) • 

Proposition 7.2. (i) The measure J is symmetric. 

(ii) For each compact subset K C T wr and open subset U such that 
Kcf/C T wr ; 



J(K x (T wr \ U)) < oo. 



22 STEVEN N. EVANS AND TYE LIDMAN 

(iii) The function Agh w1, is square-integrable with respect to J, that is, 



J J(dT',dT") A| H wr (T',T") < oo. 



Proof, (i) This is immediate from Proposition 16.31 

(ii) By construction, the measure J has the following description. Firstly, 
a weighted rooted compact M-tree T' G T wr is chosen according to P. A 
point v G T' is chosen according to the length measure px' and another 
point w G T' is chosen according to the renormalization of the weight vt> 
outside of the subtree S T ' v of points "above" v (that is, of points x such 
that v belongs to the segment [pt',x[). The subtree S T ,v is then pruned off 
and re-attached at w to form a new M-tree T" . More formally, the M-tree 
T" can be identified as the set T' equipped with new metric dy« given by 

'd(x,y), x,yeS T '' v , 

d(x,y), x,yGT'\5 T >, 

d{x, v) + d(w, y), x G S T '> V , yeT'\ S r > v , 

d(y,v) +d(w,x), y G S r > v ,x G T'\ S T '' V . 

With this identification, px" = Pt' and z^r" = V T' ■ 
We claim that if, for some e > 0, 

max dx'{v, x) < e 

x€S T '>v 

and 

v T >{S T 'n<e, 

then A GH wr (T',T") < e. Firstly, the map f :T' —* T" defined by 

( x, xeT'\S T '' v , 



d T „(x,y) :-- 



fix) 



u, x 



G S T '' 



is such that f(pT') = Pt" and 

sup{|d T /(x,y) - dr"(f(x),f(y))\ : x,y G T'} < e. 
Moreover, it is immediate that 

dp{f*VT',v T „) < e, 
and so / G F^,, T „. Note also that S T "' W = S T '' V as sets, 
max dT"(w,x) = max dT'(v,x), 

and 

v T „(s T "n = VTi(s T 'n, 

and so a similar argument shows that the map g :T" — > T' defined by 

r x, xeT"\S T '>, 

x G S r '> 
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belongs to Fj,,, T ,. Thus, Aq^i(T' ,T") < e as required. 

Now let K and U be as in the statement of part (ii). The result is trivial 
if K = 0, so we assume that K ^ 0. Since T wr \ U and K are disjoint closed 
sets and K is compact, we have that 

c:= inf A GH wr (T',T") > 0. 

T'&K,T"eU 

From what we have just observed, if < e < c, then, by Proposition 16.41 
J (if x (T wr \ [/)) < J{(T',T") : A GH wr(T', T") > e} 

< / P(dT') / fjL T '(dv)l( max d T /(u,x)>£ 
7 7t' \x£S t '>v 



+ y P(dT') ^(dv) 1 (i/(S T> ) > e 
[F + (df)2 f da 1 (max(e^) > e/2 

«60(/,a) 
+ 1 P + (d/)2y°°da ^ l(c 



< oo, 



as required. 

(iii) By an argument similar to that in part (ii) and Proposition 16.41 we have 

J j(dT',dr") A GH w r (T',r") 

= / de 2eJ {(T', T") : A G H wr (T', T") > 
7 o 

< f de2e [ P(dT') / /Wdw) 1 ( max d T '(v,x)>£ 
7o 7 7t' \x&s t ''V 



+ y de2eJv{dT') J fi T '(dv) 1 (v(S T '> v ) > e 

poo p poo , 

\ dele I P+(d/)2 / da V 1 max(e / '") > e/2 
7o 7 7o r -„ > ^ 

rco r rco 

+ de2e / P+(d/)2 / da V l(C(e^) > e) 



v€G{f,a) 

2 



= /*P+(d/)8 /" da ^ (max(e / ' u 

/poo 
P + (d/) 2 y o da 53 (C(e^)) 



< oo, 
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as required. □ 
Definition 7.3. Define a bilinear form 

£(f,g) := J J(dT',dT") (f(T") - f{T>)) (g(T") - g(T>)) , 

for /, g in the domain 

V*(£) : = {/ G L 2 (T wr ,P) : / is Borel and £(/,/) < oo}. 
Here, as usual, L 2 (T wr ,P) is equipped with the inner product 

(f,g) P := J P(dx)f(x)g(x). 

Definition 7.4. Let C denote the collection of functions / : T wr — ► M such 
that 

sup \f(T)\ <oo 

and 

\f(T')-f(T")\ 
SUP "A <T> T«\ < °° ' 

Part (i) (respectively, parts (ii) and (hi)) of the following result may be 
proved in the same manner as Lemma 7.1 (respectively, Lemma 7.2 and 
Theorem 7.3) of [EW06] . with Proposition 17.21 above playing the role played 
in [EW06J by Lemma 6.2 of that paper. Our setting is slightly different, in 
that we are working with weighted rooted compact M-trees rather than just 
weighted compact M-trees, but this doesn't require any significant changes 
in the arguments. In particular, the analogues of Lemmas 7.5, 7.6 and 
7.7 of [EW06] go through quite straightforwardly to establish the tightness 
property required for part (hi) to hold. 

Theorem 7.5. (i) The form (£,£>*(£)) is Dirichlet (that is, it is sym- 
metric, non-negative definite, Markovian, and closed). 
(ii) The set C is a vector lattice and an algebra, and C C V*{£). Hence, 
ifV{£) denotes the closure of C inV*{£)), then (£,V(£)) is also a 
Dirichlet form. 

(hi) There is a recurrent P '-symmetric Hunt process X = (X^,P ) on 
T wr with Dirichlet form {8, V {£,)). 
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